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We show that in apphcations of variational theory to quantum field theory it is essential to account 
for the correct Wegner exponent lu governing the approach to the strong-coupling, or scaling limit. 
Otherwise the procedure either does not converge at all or to the wrong limit. This invalidates all 
papers applying the so-called 5-expansion to quantum field theory. 



(T) . I. INTRODUCTION 

o . 

Variational perturbation theory is a powerful tool for extracting non-perturbative strong-coupling results from 
weak-couplingexpansions. It was initially invented in quantum mechanics as a re-expansion of the perturbation series 
. of the action m 
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which arises from splitting the potential into a quadratic part V^"^ = Q'^x'^ /2, with an arbitrary trial frequency Q, 
^ . and an interacting part 
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The perturbation expansion is then performed in powers of S, setting (5 = 1 at the end, and optimizing the result 
in 17 guided by the principle of minimal sensitivity 2J . The history and convergence properties are discussed in 
, the textbook Q. Due to the prefactor 5 in (|1.2|l . the procedure is often called (S-expansion [l(. For the anharmonic 
1-^ ■ oscillator, convergence was proved to be exponentially fast for finite J| as well as for infinite coupling strength [^IHI^- 
I ' In recent years, the method has been extended in a simple but essential way to allow for the resummation of 
divergent perturbation expansions in quantum field-theories 0, |^. The most important new feature of this field- 
theoretic variational perturbation theory is that it accounts for the anomalous power approach to the strong-coupling 
limit which the i5-expansion cannot do. This approach is governed by an irrational critical exponent lu as was first 
shown by Wegner in the context of critical phenomena. In contrast to the (S-expansion, the field-theoretic variational 
perturbation expansions cannot be derived from the action by adding and subtracting a harmonic term as in (|1.2|l . 
; ^ ' The new theory has led to the so-far most accurate determination of critical exponents via quantum field theory, as 
] amply demonstrated in the textbook In particular, the theory has perfectly explained the experimentally best 
known critical exponent a of the specific heat of the A-transition measured in a satellite orbiting around the earth 

El- ... . 

In spite of the existence of this reliable quantum-field-theoretic variational perturbation theory, the literature keeps 

offering applications of the above quantum-mechanical (5-expansion to quantum field theory, for instance in recent 
papers by Braaten and Radescu (BR) (12„ JJ5J and Ramos J^] (see also ^15j ) . 

It is the purpose of this paper to show what goes wrong with such unjustified applications, and how the proper 
quantum field-theoretic variational perturbation theory corrects the mistakes. 

II. REVIEW OF THE METHOD 

Suppose, the function /(<?) is given by a divergent series expansion around the point g = 0: 

L 

fL{9) = Y,aig\ (2.1) {WEAK} 

/=0 



2 



typically with factorial growth of the coefficients a;. Suppose furthermore, that the expected leading behavior of /(g) 
for large g has the general power structure: 
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/Af(<?) = 5" E ^"5^""' (2-2) 

m=0 

where lo is the Wegner exponent of approach to the strong-coupling limit. In quantum mechanics, this exponent 
is easily found from the naive scaling properties of the action. In quantum field theory, however, it is an initially 
unknown number which has to be determined from the above weak-coupling expansion by a procedure to be called 
dynamical determination of uj. 

Assuming for a moment that this has been done, the Lth order approximation to the leading coefficient 6o is given 
by 0: 

5(-)(.) = z-"f a,z'r^-^V^7")/-V (2.3) {BO} 
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where the z = g/fl^^°' is the variational parameter to be optimized for minimal sensitivity on z. A short reminder 
of the derivation of this formula is given in Appendix A. An application to a simple known function is shown in 
Appendix B. For a successful application to the quantum-mechanical anharmonic oscillators, the reader is referred 
to the textbook The exponent uj is equal to 2/3 for an a;^-anharmonic oscillator, and the exponentially fast 
convergence has a last term decreasing like e"'^""®'^-^^ ^ . For the oscillator, the number ui is found directly from the 
dimensional analysis in Appendix A. As mentioned above, such an analysis will not be applicable in quantum field 
theory, where to will be anomalous and must be determined dynamically. 

Most often we want to calculate a quantity f{g) which goes to a constant in the strong-coupling limit f{g) — > /*. 
This is the case for all critical exponents. Then we must set a = in H2.2|) and H2.3|) . which implies that for infinite g: 



d log g 



= . (2.4) {betaV} 



If P{g) is reexpressed as a function of /, this implies f3{f*) = 0, the standard requirement for the existence of a critical 
point in quantum field theory if f(g) = gnig) is the renormalized coupling strength as a function of the bare coupling 
strength g. 

The dynamical determination of uj proceeds now by treating not only f{g) but also the beta function H2.4(l according 
to the rules of variational perturbation theory, and determining uj to make (3* = /3(oo) vanish, which is done by 
optimizing the equation of z 



uY^-/+yi=0, (2.5) {LOG} 
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where /3i are the coefficients of the expansion of (|2.4I) with respect to g. Minimal sensitivity is reached for a vanishing 
derivative with respect to z: 



1=1 

so that z and uj are to be found as simultaneous solutions of (|2.5() and H2.6() . 



^'-i'^'j^l'/^)-^^ (2-6) {LOGl} 



III. ANOMALOUS DIMENSIONS 



As mentioned above, a number of authors have applied the (5-expansion to field theories. Most recently, this was 
done for the purpose of calculating the shift of the critical temperature in a Bose-Einstein condensate caused by a small 
interaction |l3l Il4| . Since the perturbation expansion for this quantity is a function of g/fi where /i is the chemical 
potential which goes to zero at the critical point, we are faced with a typical strong-coupling problem of critical 
phenomena. In order to justify the application of the 5-expansion to this problem, BR studied the convergence 
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properties of the method by applying it to a certain amphtude A(g) of an 0(iV)-symmetric (/)^-field theory in the 
Umit of large N, where the model is exactly solvable. 

Their procedure can be criticized in two ways. First, the amplitude A(g) they considered is not a good candidate 
for a resummation by a (5-expansion since it docs not possess the characteristic strong-coupling power structure l|2.2|) 
of quantum mechanics and field theory, which the final resummed expression will always have. The power structure is 
disturbed by additional logarithmic terms. Second, the i5-expansion is equivalent to choosing, on dimensional grounds, 
the exponent a; = 2 in (|2.2|l . which is far from the an approximate optimal value 0.843 to be derived below. Thus the 
(5-expansion is inapplicable, and this explains the problems into which BR run in their resummation attempt. Most 
importantly, they do not find a well-shaped plateau of the variational expressions A^^\g, z) as a function of z which 
would be necessary for invoking the principle of minimal sensitivity. Instead, they observe that the zeros of the first 
derivatives dzA^-'"^{g, z) run away far into the complex plain. Choosing the complex solutions to determine their final 
resummed value misses the correct one by 3% up to the 35th order. 

One may improve the situation by trying out various different w-values and choosing the best of them yielding an 
acceptable plateau in A{g,z). This happens for uj « 0.843. However, even for this optimal value, the resummation 
result never converges to the correct limit. For A{g) the error happens to be numerically small, only 0.1%, but it will 
be uncontrolled in physical problems where the result is unknown. 

Let us explain these points in more detail. BR consider the weak-coupling series with the reexpansion parameter S: 
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where a; = 



K{x)f(x)dx, 
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The geometric series in H3.1(l can be summed exactly, and the result may formally be reexpanded into a strong-coupling 
series in /i = \/l — S/{S g): 



A(5,5) 
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where 
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K{x)f-"\x)dx. 



(3.3) {BRAA2} 



The strong-coupling limit is found for h ^ where A — > fog = Jo°d,xK(x) = 1. The approach to this limit is, 
however, not given by a strong-coupling expansion of the form (|3.3|l . This would only happen if all the integrals 6^ 
were to exist which, unfortunately, is not the case since all integrals for bm with m > diverge at the upper limit, 
where 



f{x) = — arctan — ~ — . 
X 2 X 



(3.4) 



The exact behavior of A in the strong-coupling limit h —^ is found by studying the effect of the asymptotic n/x- 
contribution of f{x) to the integral in (|3.3|l . For f{x) = tt/x we obtain 



1 + h/f{x) 



7T^ + 2Tr'h-TT'h^ 
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(3.5) {BRAA3} 



The logarithm of h shows a mismatch with (|2.2() and prevents the expansion (|3.1() to be a candidate for variational 
perturbation theory. 

We now explain the second criticism. Suppose we ignore the just-demonstrated fundamental obstacle and follow 
the rules of the (5-expansion, defining the Lth order approximant A((5, cx)) by expanding H3.1|) in powers of S up to 
order S^, setting 5=1, and defining z = g. Then we obtain the Lth variational expression for 
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i(L - I + l/uj 
L-l 



(3.6) 



with w = 2, to be optimized in z. This w- value would only be adequate if the approach to the strong-coupling limit 
behaved like A + B/h' + . . . , rather than (|3.5() . This is the reason why BR find no real regime of minimal sensitivity 
on z. 

Let us attempt to improve the situation by determining uj dynamically from equation (|2.4|) . The result is w « 0.843, 
quite far from the naive value 2. This value can also be estimated by inspecting plots of A*^^)(w, /i) versus h for 
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FIG. 1: Plot of 1 — bg^'(a;,2;) versus 2 for L = 10 and uj = 0.6, 0.843, 1, 2 . The curve with cj — 0.6 shows oscihations. They 
decrease with increasing lj and becomes fiat at about ui = 0.843. Further increase of lu tilts the plateau and shows no regime 
of minimal sensitivity. At the same time, the minimum of the curve rises rapidly above the correct value of 1 — 6o = 0, £is can 
be seen from the upper two curves for a; = 1 and w = 2, respectively. 




0.5 1 1.5 12 




0.5 1 1.5 



FIG. 2: Left-hand column shows plots of 1 - b}, '{uj,z) for L = 10, 17, 24, 31, 38, 45 with w = 2 of <5-expansion of BR, 
right-hand column with optimal uj = 0.843. The lower row enlarges the interesting plateau regions of the plots above. Only 
the right-hand side shows minimal sensitivity, and the associated plateau lies closer to the correct value 1 — 6o = than the 
minima in the left column by two orders of magnitude. Still the right-hand curves do not approach the exact limit for L — » oo 
due to the wrong strong-coupling behavior of the initial function. 
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FIG. 3: Deviation of 1 
reached, missing the correct number by about 0.1%. 



^o^piatcau ("^ = 0.843) from zero as a function of the order L. Asymptotically the value —.001136 is 



several different w-values in Fig. ^ and selecting the one producing minimal sensitivity. It produces reasonable 
results also in higher orders, as is seen in Fig. [3 The approximations appear to converge rapidly. But the limit 
does not coincide with the known exact value, although it happens to lie numerically quite close. Extrapolating 
the successive approximations by an extremely accurate fit to the analytically known large-order behavior p] with a 

function &o^piatcau ('^ ~ 0-843) = A + B L~'^, we find convergence to ^ = 1 — 0.001136, which misses the correct limit 
A=l. The other two parameters are fitted best by S = -0.002495 and k = 0.922347 (see Fig.|g. 

We may easily convince ourselves by numerical analysis that the error in the limiting value is indeed linked to 
the failure of the strong-coupling behavior (|3.5|l to have the power structure H2.2|) . For this purpose we change the 
function f{x) in equation (|3.2|1 slightly into f{x) f{x) — f{x) + 1, which makes the integrals for 5,„ in 1)3. 3|) 

convergent. The exact limiting value 1 of A remaines unchanged, but bQ^"* acquires now the correct strong-coupling 
power structure (|2.2I) . For this reason, we can easily verify that the application of variational theory with a dynamical 
determination of u yields the correct strong-coupling limit 1 with the exponentially fast convergence of the successive 
approximations for L — > oo like 6q'^'' w 1 — exp (—1.909 — 1.168 L). 

In the next section we are going to point out, that an escape to complex zeros which BR propose to remedy the 
problems of the (5-expansion is really of no help. 



IV. THE MYTH OF COMPLEX ZEROS AND FAMILIES 

It has been claimed |l6j| and repeatedly quoted that the study of the anharmonic oscillator in quantum 
mechanics suggests the use of complex extrema to optimize the (5-expansion. In particular, the use of so-called 
families of optimal candidates for the variational parameter z has been suggested. We are now going to show, that 
following these suggestions one obtains bad resummation results for the anharmonic oscillator. Thus we expect such 
procedures to lead to even worse results in field-theoretic applications. 

In quantum mechanical applications there are no anomalous dimensions in the strong-coupling behavior of the 
energy eigenvalues. The growth parameters a and lu can be directly read off from the Schrodinger equation; they 
are a = 1/3 and lu ^ 2/3 for the anharmonic oscillator (see Appendix A). The variational perturbation theory is 
applicable for all couplings strengths g as long as h\^\z) becomes stationary for a certain value of z. For higher orders 
L it must exhibit a well-developed plateau. Within the range of the plateau, various derivatives of b\^\z) with respect 
to z will vanish. In addition there will be complex zeros with small imaginary parts clustering around the plateau. 
They are, however, of limited use for designing an automatized computer program for localizing the position of the 
plateau. The study of several examples shows that plotting h\^\z) for various values of a and lo and judging visually 
the plateau is by far the safest method, showing immediately which values of a and oj lead to a well-shaped plateau. 

Let us review briefly the properties of the results obtained from real and complex zeros of dzh\^\z) for the anhar- 
monic oscillator. In Fig. ^ the logarithmic error of ^g^'' is plotted versus the order L. At each order, all zeros of the 
first derivative are exploited. To test the rule suggested in , only the real parts of the complex roots have been 
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used to evaluate 6q^^ . The fat points represent the results of real zeros, the thin points stem from the real parts of 
complex zeros. It is readily seen that the real zeros give the better result. Only by chance may a complex zero yield 
a smaller error. Unfortunately, there is no rule to detect these accidental events. Most complex zeros produce large 
errors. 

We observe the existence of families described in detail in the textbook 3] and rediscovered in Ref. These 
families start at about = 6, 15, 30, 53, respectively. But each family fails to converge to the correct result. Only 
a sequence of selected members in each family leads to an exponential convergence. Consecutive families alternate 
around the correct result, as can be seen more clearly in a plot of the deviations of 6q^'' from their L —> oo -limit 

in Fig. [SI where values derived from the zeros of the second derivative of b^^^ have been included. These give rise 
to accompanying families of similar behavior, deviating with the same sign pattern from the exact result, but lying 
closer to the correct result by about 30%. 
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FIG. 4: Logarithmic error of the leading strong-coupling coefficient of the ground state energy of the anharmonic oscillator 
with a;* potential. The errors are plotted over the order L of the variational perturbation expansion. At each order, all zeros 
of the first derivative have been exploited. Only the real parts of the complex roots have been used to evaluate 6g^'. The fat 
points show results from real zeros, the smaller points those from complex zeros, sizedecreasing with distance from real axis. 



V. TEMPERATURE SHIFT FOR N = 2 REVISITED 

Much attention has been paid to a field theoretic model with 0(2)-symmetry to calculate in a realistic 

context the coefficient ci, which enters into the temperature shift of the Bose-Einstein condensation parametrized as: 

AT, 



(0) 



Cian 



1/3 



(5.1) {Neq2A} 



Presently, five coefficients of the relevant perturbation expansion are known for the weak-coupling expansion [l3LIT3 . ll^ 

3 

F(x) = ^ a„ , (5.2) {Neq2B} 

n=-l 

whose asymptotic value for x ^ oo coincides with ci: ci — F* = lim^^oo F{x) ■ The known coeflicients are a_i = 
-13.9707, ao = 0, ai = -0.446572, az = 0.264412, = -0.199. 

We would like to offer an alternative resummation result for this series to that in Ref. ^3 • It is based on considering 
the function x F{x) containing no negative powers of x. The desired number ci is the the leading coefficient bo of the 
strong-coupling expansion 



F{x) = X 



(5.3) {Neq2D} 
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FIG. 5: Deviation of the coefRcient 6q^' from the exact value is shown as a function of perturbative order L on a linear scale. 
As before, fat dots represent real zeros. In addition to Fig. |1] the results obtained from zeros of the second derivative of 6q^' are 
shown. They give rise to own families with smaller errors by about 30%. At A'' = 6, the upper left plot shows the start of two 
families belonging to the first and second derivative of feQ^', respectively. The deviations of both families are negative. On the 
upper right-hand figure, an enlargement visualizes the next two families starting at A = 15. Their deviations are positive. The 
bottom row shows two more enlargements of families starting at A = 30 and TV = 53, respectively. The deviations alternate 
again in sign. 



The result for ci should be unaffected by this modification of the function, and given by the optimized ith-order 
approximations 

4^)(z,C.) = X:«.-'-^('^"' + ^_':'^/"). (5.4) {Neq2E} 

1=0 ^ ^ 

For the available orders £ < 4, this set of functions is now inspected for plateaus. For L < 3 there is none. For L = 3 
and L — A, a, plateau can be identified unambiguously as the only horizontal turning point solving simultaneously 
dzc[^\z,uj) = and d'^c[^\z,uj) = 0. The results are 

L = 3 z(^) = 1.089 w(^) ^ 1.071 cf^ = 0.940 (5.5) {cl} 

L = A z^^^ = 2.057 w^"*) 0.571 c^^^ ^ 1.282 (5.6) 

Given only two approximations for Ci it is unrealistic to attempt an extrapolation to L — > oo, as done with another 
selection rule of optima in Ref . , but it is interesting to note that the value of the coefficient for the temperature 
shift c^^'^ = 1.282 is in excellent agreement with the latest Monte Carlo result of ci « 1.30 [l9|. 
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VI. RENORMALIZATION GROUP AND VARIATIONAL PERTURBATION 

The most convincing evidence for the power of the field-theoretic variatio nal p erturbation theory with anomalous 
dimensions comes from applications to critical exponents in 4— e dimensions [^llOj. The results obtained turn out to be 
immediately resummed expressions of the e-expansions, which can be recovered as a Taylor series. The renormalization 
group function f3{g) is obtained from the weak-coupling expansion of the renormalized coupling constant g in terms 
of the bare coupling constant gs 0, 113 : 

»t \ d\ogg{gB,e) rd log 35(3, e)] -i 

f3{g,ej = -e g —r-. = "£5 — ■ (6-1) {beta} 

d log gs L d log g J 

Due to renormalizability, P{g) necessarily has the form 

(3{g,e)^~eg + f^o{g) . (6.2) {beta2} 

Perturbation theory with minimal subtractions yields the weak-coupling expansion: 



9 = 9B + J2 e , (6.3) {g-Coupl: 



k=l 



where fkigs) possesses an expansion in powers oi gs, starting with gs^^^ ■ By suitably normalizing g and gs, the 
leading coefficient of /i can be made equal to minus one: fiigs) = —gs^ + 0(173)'^. The function /?o(5) can be 
expressed in terms of the residue fiigs) of the e-pole in eauation (|6.3|) alone: 

f3o{g) = fi{g) - g f[{g) ■ (6.4) {betao} 

Recall the standard proof for this based on combining Eqs. Htj.lfl and Ht).2|l to 

Po{g) ^ eg ~ egs g{gB,e) , (6.5) {bx} 

which becomes, after inserting equation (|6.3(l : 

00 00 

0o[gB + e-''] - e ^ [Mgs) - .9^/^55)] e"' • (6.6) {bxx} 

fc=l k=l 

The limit e — > cx) leads directly to the property Ht).4|l . 

Another well-known fact is that all the functions fkigs) for fc > 1 can be expressed in terms of the residues fiigs) 
only 10]. Indeed, taking the derivatives of /3o(5) in equation (|6.5|) with respect to gs and e: 

f3'o{g)dg,g^-egBdlg, (6.7) {by} 

P'oia) deg ^ g + ed^g- gs dg^g -egs dg^d^g , (6.8) 

eliminating (^^{g) between these two equations, and inserting the expansion (|6.3|l . we obtain order by order in 1/e a 
recursive set of differential equations for the functions fkigs) with k > 1, which are power series in gs- If we now 
expand 

00 00 

higs) = -gs'^ + ^ijgB' , fkigs) = '^^jgn^, (6.9) {power} 

J=3 j=k+l 



a solution is readily found, beginning with 



7fe,fc+i = 


(-1)', 




12A = 


73,5 = 


29 

y73, 




72,6 = 


74,6 = 


37 

-y73, 




72,7 = 


74,7 = 


751 2 


141 


75,7 = 



8 3 7 

-373, 72,5 = 2 73 - 2'>'4' (6-10) {soli} 

18 22 32 2 39 ,^ , 

y 73 74 -y 75, 73,6 = — ^73 + y74, (6.11) 

n 2 13 16 5 3 551 59 

274 + y 73 75 - y 76, 73,7 = 2 ^3 - y^- 7374 + y 75, (6.12) 

103 

— 73. (6.13) {S0I2} 
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In the renormalization group approach, a fixed point ^ is determined by the zero of the /3-function: (3{g*) = 0. 
The Wegner exponent uj governing the approach to scaUng is given by the slope at the fixed point: oj — (3'{g*). The 
two quantities have e-expansions 
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The coefficients aj and lOj are determined from the residues 7j as: 
ai = 1, a2 = 2 73, 



"3 = 8 7^ + 3 74 



a4 = 40 7f + 30 73 74 + 4 75, 



as = 224 7* + 252 7^ 74 + 27 7I + 48 73 75 + 5 76, 



and 



UJ2 



'2 73, 



W3 



-8 7^ - 6 74, 



-40 7I - 48 73 74 - 12 75- 



(6.14) {gStar} 

(6.15) {alpha} 

(6.16) {alpha2} 

(6.17) {omega} 



We can now convince ourselves that precisely the same results can be derived from variational perturbation theory 
applied to the weak-coupling expansion (|6.3I) (and as shown in l2l| from the expansion of any other critical exponent). 
We determine uj dynamically solving Eq. H2.5|l . We insert for lo an unknown e-expansion of the form H6.14|l . The 
variational parameter z is then adjusted to make (|2.5|l stationary. Then, since for e — > the weak-coupling coefficients 
of gigs) in the expansion (|6.3|) behave like ~ e^~', z has to scale with e, so that we may put z — Cie+C2 £'^+C3^^+0{e*), 
and solve equations H2.5() and (|2.6|l for each perturbative order L, order by order in e. This leads to a rapidly increasing 
number of non- linear and not even independent equations for the unknown Q and lui , some depending also on the 
order L. 

Despite these possible complications, the solutions turn out to be well structured and easily obtained. At each L 
to lowest order in e, the term independent of e in H2.5|l and the coefficient of in H2.6|l demand that Ci = 1- In 
addition, they require jk.k+i = ("1)*^ for some k, in agreement with Eqs. (|6.10|) . Such conditions imposed on jk.i 
can, of course, not depend on the order L, but must be enforced in general. Raising the order of e in (|2.5|l and (|2.6|l . 
and imposing ^1 = 1 as well as the conditions already established for the "fk,ii all dependences on the ujk and Cfe 
disappear, and we are left with conditions on 7^^; alone, which reproduce exactly the relations (|6.10() through H6.13|l . 
This shows, that the variational perturbation method is completely compatible the well-known e-expansions, if the 
input divergent series has a structure satisfying the renormalization group equation (|6.2|l . 

After having reproduced "fk,ii there are further equations to be solved. Going to the next higher order in e, cither 
for H2.5|l or for (|2.6I) . gives a relation involving exactly one of the expansion coefficients of e/w, which are simply 
related to the coefficients uji of uj. In this way, the renormalization group results of H6.17|l are exactly reproduced. 
These solutions are stable in the sense, that with increasing order L, the expansion coefficients ui for I < L remain 
unchanged. This proves, that the variational method produces the same e-expansions of all critical exponents as 
renormalization group theory. At the same time this implies that the standard (5-expansion which does not allow for 
the anomalous dimension w is bound to fail. 

It is noteworthy, that several other conditions are automatically satisfied up to some order e^, e^~^, or e^"^, 
respectively. Among them is the variationally transcribed second logarithmic derivative of the weak-coupling series 
and the derivative thereof: 



E ifL-l + l/uj 
h,z^ \=-l-uj, 

(=1 



where the hi are the expansion coefficients of 



9B 9" {93) 
9'(.9b) 



(6.18) 
(6.19) 

(6.20) 



Of some computational benefit is the observation, that with the same accuracy in e the first and second derivatives of 
the variational series H2.3|l themselves vanish (here for a = 0). This means, that the function has a flat plateau. For 
a typical field-theoretic application with only a few known perturbation coefficients, the plateau is easily found by 
inspection. Therefore, if the model possesses a well-behaved /^-function satisfying equation (|6.2|l . we expect a reliable 
result for the anomalous dimension w if it is chosen such as to produce an acceptable plateau. The ordinate of the 
plateau is the most promising variational perturbative value for the quantity analyzed to the respective order. 
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VII. CONCLUSION 

Summarizing this paper we have learned that the so-caUed (5-expansion is inapphcable to quantum field theory, 
since it does not account for the Wegner exponent oj of approach to the strong-coupling limit. Only the field-theoretic 
variational perturbation theory yields correct results by incorporating w in an essential way. 

VIII. APPENDIX A 

Here we review briefly how the strong-coupling parameters a and w in 1)2. 2|l and the variational equation l|2.3() 
for the leading strong-coupling coefficient are found for the anharmonic oscillator with the Schrodinger equation in 
natural units 

-i*" + y«' + .ga;2«^' = £;*. (8.1) {Al} 

We rescale the space coordinate x so that the potential becomes 

V{x) = ig-'/("+i^a;2 + , (8.2) {A2} 

any eigenvalue has the obvious strong-coupling expansion 



iJ = gV(^+i)^6,(g-V(^+i))', (8.3) {A3} 

1=0 

where hi are the strong-coupling coefficients. The aim is to determine them from the known weak-coupling coefficients 
an of the divergent perturbation expansion: 

oo 

E^Y.^'ia'- (8.4) {A4} 

1=0 

The solution of this problem comes from physical intuition, suggesting that the perturbation expansion should be 
performed around an effective harmonic potential r2^a;^/2, whose frequency is different from the bare value 1/2 in 
()8.1|l . depending on g and the order L of truncation of H8.3|l . Thereafter only the difference between the anharmonic 
part and the effective harmonic part is to be treated by perturbation methods. The trial frequency of the effective 
potential can be fixed later by the consideration, that the resulting quantity of interest should be as independent as 
possible of SI, according to the principle of minimal sensitivity. With the harmonic trial potential Vq^^ = r2^a;^/2, the 
interaction potential H1.2|l reads Vq^ — 5 \gx'^'^ — (Sl^ — l)x^/2]. The parameter 5 organizes the reexpansion and is 
set equal to 1 at the end. The expansion proceeds from the rescaled Schroedinger equation H8.1() : 

I „ x"^ Sqx'^'^-^ E , , 

^2*" + y* + ^^ = ^*, (8.5) {A5} 



where fi = -^/Sl^ — (5(S1^ — 1). To order L, the energy has the reexpansion 

E^'^K^.9)-Pj:af\-1^)\ (8.6) {A6} 

with the well-known weak-coupling expansion coefRcients as defined in equation (|8.4(l . The strong-coupling behavior 
()8.3(l suggests changing the variational parameter from SI to z := jprrr- In the limit g ^ oo we obtain the reexpansion 
which must be optimized in z: 

E^-\z)=9-j:af^z^-(''-'\^'-''^'^) (8.7) {A8} 

1=0 ^ ^ 



where w = 2/(k + 1) and a = 1/(k -|- 1). For the leading coefficient of the strong coupling expansion of the ground 
state energy, Ea. H8.7|l leads directly to the variational equation (|2.3|) . 
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IX. APPENDIX B 



In order to gain further insight into the working of the variational resummation procedure, we apply it to the simple 
test function 

/(a;) = (l + x)"=x"(l + i)" (9.1) 

with weak coupling coefficients a„ — (^) and a leading strong-coupling behavior / ^ x°'{l + a/x + . . .), so that &o = 1- 
Inserting this information into equation H2.3fl . we obtain the variational leading coefficient to Lth order: 

'rw=|:(:)(r_:)='-". 

which is easily transformed into the expression: 



Determining the variational parameter z according to the principle of minimal sensitivity requires a well developed 
plateau of &q^'' as a function of z. For the simple test function, the derivative dzh\^\z) can be obtained in the closed 
form: 

L 



This exhibits a flat plateau around z = 1 if the order L is much larger than a. An equally flat plateau is found for 
6q^'(z). The value of the leading strong coupling coefficient 6q^' at the plateau is 

C(i)^f°)i:(^)^(-i)-^i, («) 

(=0 



in perfect agreement with the exact result, thus confirming the applicability of the resummation scheme for this class 
of problems. 
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